PROSPECTS FOR HELIO- AND ASTEROSEISMOLOGY 



J0rgen Christensen-Dalsgaard 

Institutfor Fysik og Astronomi, Aarhus Universitet, DK-8000 Aarhus C, Denmark 



ABSTRACT 



Major progress has been made in helio- and asteroseis- 
mology in recent years. In helioseismology, much of the 
activity has been in local helioseismology. However, the 
recent revision of solar surface abundances, and the re- 
sulting problems in reconciling solar models with the he- 
lioseismic inferences, have lead to renewed activity in so- 
lar modelling and global helioseismology. Interesting, al- 
though perhaps not compelling, evidence has been found 
for solar g modes in observations with the GOLF instru- 
ment on the SOHO spacecraft. Extensive asteroseismic 
results have been obtained from ground-based observa- 
tions as well as from the WIRE and MOST satellites, 
and much is expected from the upcoming launch of the 
CoRoT satellite and, in a few years, from the Kepler 
mission. In parallel, stellar modelling is being extended 
to take some account of hydrodynamical effects, while 
large-scale hydrodynamical calculations are providing in- 
creasingly realistic simulations of these effects. The out- 
come of these activities will undoubtedly be a far bet- 
ter understanding of stellar internal properties and stel- 
lar evolution, together with an improved insight into the 
physics of matter under the extreme conditions found in 
stars. 
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1. INTRODUCTION 



The conference demonstrated the impressive and exciting 
progress that has been made over recent years in the field 
of helio- and asteroseismology, as well as the challenges 
that we face and the promising prospects of overcom- 
ing them. Global helioseismology is a mature field, and 
progress is unavoidably slower than in the early stages. 
The continuing accumulation of data is essential to study 
the intriguing variations with time that have been found 
in the Sun's internal dynamics but will only very slowly 
lead to further significant improvements in the data on the 
Sun's average properties. Such improvements are likely 
to come rather from improved techniques for the analysis 
of the helioseismic data. At the same time we are faced 



with increasing challenges to our understanding of solar 
modelling. Local helioseismology is still under rapid de- 
velopment, with an encouraging increasing emphasis on 
the understanding of the diagnostic potential and relative 
merits of the different techniques. Asteroseismology has 
a long and fruitful history in the study of some classes 
of stars, particularly the compact pulsators. However, the 
development of instruments and the organization of ma- 
jor campaigns over the last few years have led to an explo- 
sion of activity on the observational front, with the efforts 
towards modelling and data interpretation attempting to 
keep up. In all these areas dramatic progress is expected 
from new satellite missions, and other instrumentation, to 
be launched over the coming few years or in the planning 
stage. 

In the present paper I admit to a strong bias towards 
global helioseismology and asteroseismology (although 
not entirely to the 'dark side', as Jesper Schou once ac- 
cused). Indeed, it is obvious that asteroseismology is a 
natural extension of helioseismology in the attempt to un- 
derstand the large-scale structure and dynamics of stel- 
lar interiors and their evolution with time; the restriction 
to far fewer modes in asteroseismic investigations is to 
some extent compensated by the ability to observe stars 
over a broad range of parameters. To complement this 
biased view, Toomre (these proceedings), in his intro- 
duction to the conference, provides an overview of the 
accomplishments and promise of local helioseismology 
in the study of the more complex three-dimensional and 
strongly time-dependent phenomena in the outer parts 
of the Sun, including the aptly named 'solar sub-surface 
weather' . Also, Birch (these proceedings) provides an ex- 
cellent review of the techniques of local helioseismology. 

In the analysis of the helio- and asteroseismic data we 
should keep in mind that the ultimate scientific goal is to 
obtain insight into and understanding of the many com- 
plex phenomena that we see in stars. As tools in this pro- 
cess, and to emphasize the excitement of the field, we 
obviously use images illustrating the results, but such im- 
ages cannot be goals in themselves. (And indeed, through 
being graphically striking, they might even occasionally 
mislead our intuition, unless we have a firm understand- 
ing of how the graphical representation relates to the un- 
derlying solar or stellar properties.) It is very encourag- 
ing that the conference showed several examples of a shift 
towards a comparison of the observational data with, one 
may hope, realistic models of phenomena under investi- 



gation. This is an essential part of relating the observa- 
tions to the physics of the stellar interiors. 

The conference took place just over a year after George 
Isaak sadly passed away. George played a seminal role 
in the development of helioseismology, through the ini- 
tial observations of global five-minute oscillations and 
the establishment of the BiSON network and the BiSON 
group; he furthermore strongly encouraged the extension 
of these studies to other stars. Thus it was very appropri- 
ate that the section of the conference dealing with low- 
degree seismology from the Sun to the stars was dedi- 
cated to his memory. He would have greatly enjoyed tak- 
ing part in the conference, with 'a comment and a ques- 
tion' after many of the talks, no doubt often exposing 
weaknesses in the argument or pointing the way forward. 



2. IN FOR THE LONG HAUL 



As pointed out by Howe (these proceedings) we now have 
helioseismic data spanning more than 30 years, or 1 1/2 
solar (magnetic) cycles. The BiSON and GONG net- 
works and the instruments on the SOHO satellite have 
covered a full 11 -year sunspot cycle with data with a 
very high duty cycle. As discussed by Howe, these re- 
sults have shown intriguing variations in the properties of 
the solar interior and the physics of the modes, in most, 
but not all, cases clearly related to the sunspot cycle. 
The variation of the average multiplet frequencies, or the 
component of frequency splittings even in the azimuthal 
order m, is generally closely correlated with the surface 
magnetic field ifTTl Uj[ l6lll . However, in his presenta- 
tion Verner showed that for low-degree modes there is a 
clear hysteresis between the frequency shift and various 
indicators of solar activity (Chaplin et al. . the se proceed- 
ings), as has previously been found 164. Il04ll . The long 
time basis in the results presented by Verner allowed the 
demonstration of interesting variations in the relation be- 
tween frequency shift and activity between solar cycles. 

The solar internal rotation shows more striking variations 
with time. The clearest signal are the zonal flows, re- 
gions of slightly slower and faster rotation (Howe, these 
proceedings). At low and intermediate latitude a similar 
pattern in the solar surface velocity has previously been 
identified as solar 'torsional oscillations' 16011 : from he- 
lioseismic analyses it is now known that these flows ex- 
tend over a substantial fraction of the convection zone. 
They converge towards the equator as the sunspot cycle 
progresses and indeed the location of the flow is close 
to, but shifted slightly towards lower latitudes than, the 
sunspot belts, possibly suggesting that the flows lead the 
activity in the development of the solar cycle Il62ll . At 
higher latitudes the corresponding surface feature had 
not previously been identified; there the bands of slower 
and faster rotation move towards the pole with the pro- 
gression of the sunspot cycle. The physical causes for 
these flows are discussed by Rempel (these proceedings); 
the high-latitude component appears to be dominated by 
the feedback of the Lorentz force on the differential ro- 



tation, whereas no entirely satisfactory explanation has 
been found for the behaviour at lower latitudes. 

As discussed by Howe (these proceedings) an oscillation 
with a period of 1.3 yr has been detected near the base of 
the convection zone, at the equator. This was observed 
consistently in independent analyses during the period 
1996 - 2001 while the behaviour in recent years has been 
more erratic. Basu & Antia IU8I1 . while seeing a superfi- 
cially similar behaviour in their results, argued that it was 
not significant, and indeed the rather erratic nature of the 
observed variation might reasonably raise concerns about 
the physical reality of the 1.3-yr oscillation. On the other 
hand, it is not unreasonable that oscillations in a poten- 
tially strongly coupled and nonlinear system display time 
variations, either tied to specific phases in the solar cycle 
or in a more chaotic fashion. 

Given these time variations, strongly related to the solar 
cycle and of a somewhat irregular nature, it is obviously 
important to continue detailed helioseismic analyses on a 
timescale covering several solar cycles, and in a form that 
ensures homogeneous datasets. For the ground-based ob- 
servations this requires maintaining both the BiSON and 
GONG networks; it is important to emphasize that these 
two networks complement each other, the BiSON net- 
work providing data on the lowest-degree modes that are 
superior to those obtained by the spatially resolved ob- 
servations. A further important incentive for extending 
the BiSON observations is the push to lower frequencies 
in the p-mode range. These modes have lifetimes of sev- 
eral years and hence their frequencies can be determined 
with extremely high precision; thus they are potentially 
very valuable for the determination of the core structure 
and rotation. The space-based observations clearly have 
a limited lifetime; however the GOLF observations on 
SOHO should be continued as long as possible, and a rea- 
sonable overlap must be ensured between the MDI obser- 
vations on SOHO and the HMI observations on SDO, to 
ensure an adequate intercalibration of these two datasets. 



3. IMPROVED DATA FOR GLOBAL HELIOSEIS- 
MOLOGY 

Global helioseismology has been remarkably successful 
in uncovering the detailed structure and rotation of the 
solar interior. Yet, unavoidably, these inferences have 
revealed a need for even finer details, to allow a deeper 
understanding of the underlying phenomena. Results on 
the structure of the inner core are ambiguous, at a level 
which may have some effects on the use of solar models 
to constrain neutrino properties. The determination of the 
rotation of the core is rather uncertain, leaving room for a 
small core rotating substantially faster (or perhaps more 
slowly) than the bulk of the radiative interior. The im- 
portant issues of penetration and overshoot at the base of 
the convection zone, not least the potential variation with 
latitude, have only been addressed in fairly crude terms. 
Also, substantial uncertainty remains concerning the de- 
tailed properties of the so-called tachocline, i.e., the tran- 
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Figure 1. Model of a one-year m-averaged power spec- 
trum from the MDI instrument on the SOHO spacecraft, 
for spherical harmonic degree I = 6 and radial order 
n = 12. The thin line shows the d ata an d the heavier 
grey line shows the fit. Adapted from hlO&l . 

sition region between the latitudinally varying rotation in 
the convection zone and the near-constant rotation in the 
radiative interior, and its location relative to the base of 
the convectively unstable envelope. 

Needless to say, resolving these issues requires better 
data; in particular, not only random, but also potential 
systematic, errors must be reduced. In principle, random 
errors can be reduced through extending the observations 
i67il . although a significant improvement at the present 
stage will obviously require perhaps unrealistically long 
time series (at least on the natural timescale of the present 
author). Also, the variation of the frequencies over the so- 
lar cycle makes it difficult to combine data over such long 
periods. Indeed, the study of these time variations is in it- 
self of great interest and clearly requires a determination 
of accurate frequencies from time series of limited dura- 
tion. Thus it is of obvious interest to develop techniques 
that allow reliable frequency determinations, or perhaps 
more direct inferences about the solar interior, from the 
existing observations. 

The analysis of helioseismic observations is complicated 
by the inability, owing to the restriction of the observa- 
tions to one hemisphere of the Sun, to separate fully the 
contributions from modes of different degree I and az- 
imuthal order m in the spatial analysis. Thus the analysis 
of the time series for a given target (Iq, tuq) must take into 
account the cross-talk from modes of other I and m. In 
earlier analyses this has been done in a somewhat incom- 
plete way IBl l90t 19111 . Also, damping and the stochas- 
tic nature of the excitation of the modes must be taken 
into account, leading, for example, to an asymmetry in 
the peaks in the power spectra. If not treated properly, 
such effects may cause systematic errors in the frequen- 
cies which could lead to serious problems in the infer- 
ences of the solar internal properties. 

Substantial progress towards resolving these problems 
have been made by Jefferies, Vorontsov & Giebink (these 



proceedings), who report results from a major project to 
develop techniques for analysis of global helioseismic 
data which, as far as possible, take the known proper- 
ties of the modes into account. The cross-talk is take n 
fully into account, using techniques developed in 1 1 08 ] . 
Also, the physics of the excitation and damping is de- 
scribed in some detail, resulting in a good representation 
of the spectral line profiles. Finally, the solar structure 
is directly included in the fit, through a phase function 
corresponding to the inferred solar seismic model. An 
example of the resulting fit to the observed power spec- 
trum is illustrated in Figure ^ The results obtained by 
Jefferies et al. show a marked reduction in the frequency 
residuals resulting from the fit, compared with the previ- 
ously available frequencies, and provide very interesting 
inferences of the sound-speed correction relative to a so- 
lar model and the buoyancy frequency, even in the inner 
solar core. 

The analysis by Jefferies et al. considered mean multi- 
plet frequencies and concentrated on the structure of the 
solar interior. A major challenge will be to extend the 
analysis to individual frequencies, in this way determin- 
ing the rotational splitting and possible structural depar- 
tures from spherical symmetry. Another important is- 
sue, particularly for the study of the solar core, is the 
consistent combination of different datasets, particularly 
the spatially resolved data from GONG or MDI with the 
more detailed and accurate low-degree data from the un- 
resolved observations with BiSON or GOLF. Such tasks 
must have very high priority. There seems little doubt that 
the full development of these techniques will very consid- 
erably strengthen our ability to infer the interior structure 
and dynamics of the Sun, including the determination of, 
or constraints on, possible time variations. 

High-degree modes are particularly sensitive to instru- 
mental effects and, for ground-based observations, effects 
of the Earth's atmosphere. Thus, despite their strong po- 
tential as diagnostics of the outer parts of the convection 
zone 1 79], they have so far seen little use in global helio- 
seismology. It is encouraging that Rabello-Soares et al. 
(these proceedings) report progress on the understanding 
and elimination of instrumental effects for the MDI in- 
strument on SOHO. If combined with the analysis tech- 
niques developed by Jefferies et al. we may hope finally 
to obtain reliable data on high-degree modes for use in in- 
verse analyses. This will be particularly important for the 
study of the thermodynamics of the hydrogen and helium 
ionization zones. 



4. SOLAR G MODES? 



Solar g modes have been of the nature of a 'holy grail' 
of helioseismology ever since the early claims to the de- 
tection of a 160-minute oscillation [25, 94, 88]. The rea- 
son for this interest is obvious: g modes have their largest 
amplitudes in the solar core and hence potentially provide 
much more detailed information about the structure and 
rotation of the deep solar interior. This is perhaps most 



important in the case of rotation, where no information is 
obtained from radial modes which penetrate most deeply, 
and where the low-degree nonradial p modes have very 
limited sensitivity in the inner core. Consequently, there 
have been extensive efforts and further claims to the de- 
tection of g modes 11371 14711 - There seems little doubt that 
the early claims were premature; a stringent upper limit 
to the amplitude of possible g modes, well below the ear- 
lier apparent detections, was obtained in |6]. 

Theoretical predictions of the g-mode amplitudes are ex- 
tremely uncertain but tend to support this non-detection. 
Gough [51] and Kumar et al. |66] predicted amplitudes 
around 0.1 cms -1 and 0.01 cms" 1 , respectively, at a fre- 
quency of 200 /iHz, with a rapid decrease of amplitude at 
lower frequencies. These estimates were based on ex- 
tending the stochastic excitation by near-surface convec- 
tion, which is quite successful for the p modes (cf. Sec- 
tion!^}, to lower frequencies. Excitation of gravity waves 
also likely occur at the base of the convection zone (see 
Talon, these proceedings), although most efficiently at the 
longer periods corresponding to the characteristic con- 
vective time scales in the lower parts of the convection 
zone. In any case, these estimates are probably suffi- 
ciently uncertain that they cannot be taken as strict re- 
quirements for an observational search for g modes. 

Due to their potential diagnostic importance, g-mode de- 
tection has been an important goal of the GOLF instru- 
ment on the SOHO spacecraft [48]. The analysis of the 
GOLF data has in fact led to a renewed interest in the 
search for solar modes at low frequencies, including g 
modes. Turck-Chieze et al. found tentative evi- 

dence for multiplet structure at very low amplitudes, at 
frequencies corresponding to low-order p and g modes. 
In these proceedings, Garcia et al. present an analysis that 
instead attempts to identify the expected overall structure 
of high-order g modes. For such modes the periods are 
approximately uniformly spaced, at a spacing that de- 
pends on the degree. By detecting this structure in the 
observations the presence of g modes can be confirmed, 
and the properties of the modes, including the average pe- 
riod spacing and possibly individual frequencies, can be 
determined. Such analyses have been extensively used 
in earlier searches for, and claimed detections of, solar 
g modes; however, the analysis by Garcia et al. probably 
represents the first time that they are applied to the GOLF 
data. 

Garcia et al. analyse a 10-yr time series of GOLF data; 
they remove a solar background from the power spec- 
trum and carry out a Fourier transform of the resulting 
spectrum, as a function of period, to identify possible pe- 
riodic structures, as expected from asymptotic theory. To 
ensure that asymptotic theory is valid, while the expected 
amplitudes are not too small, the analysis is restricted to 
the cyclic frequency interval between 25 and 140 /iHz. 
Analysis of corresponding artificial data based on a solar 
model shows a strong signal only at the period spacing 
corresponding to I = 1. As illustrated in Figure [2] the 
observed frequencies show a similar response at approx- 
imately the same period spacing as for the model. 
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Figure 2. Power spectrum of power-spectrum density in 
the g-mode frequency range. The lower panel is for a 
solar model, while the upper panel is based on 10 years 
of observations with the GOLF instrument on the SOHO 
spacecraft. The shaded area indicates the location of the 
peak corresponding to dipolar g modes equally spaced in 
period. (From Garcia et al., these proceedings.) 

Given the extensive analysis required to obtain this result, 
and the potentially complex sources of noise, including 
solar 'noise', it is evidently important to test its signifi- 
cance of the result. Based on a Monte Carlo simulation 
with a large number of noise realizations Garcia et al. 
find that the observed / = 1 peak is significant at a level 
of 99.5 %. This is encouraging; however, further tests are 
probably required using different assumptions about the 
character of the noise background. As a related issue, it 
is found that the observed spectrum of the spectrum can 
be reproduced in the presence of background noise only 
if the lifetimes of the modes are relatively short, of the 
order of a few months. This is substantially below the ex- 
pected damping times of these modes; however, it could 
reflect phase variations of the modes on this timescale, 
and caused by variations in the internal properties of the 
Sun. The possible 1.3-yr oscillation in the rotation rate 
indicates that such variations may be present. 

From model fits to the observed data Garcia et al. con- 
clude that the core rotation of the Sun may exceed the 
surface value by as much as a factor five. This would ob- 
viously be a very interesting result but probably requires 
further confirmation, as well as tests of whether it is con- 
sistent with existing constraints from p-mode splittings 
on the rotation of the deep solar interior |26]. 

Unfortunately, the analysis does not provide estimates of 
the oscillation amplitudes, and hence comparisons can- 
not be made with earlier upper bounds on the g-mode 
amplitudes or with the theoretical predictions. Such es- 
timates can presumably be obtained from further simula- 
tions, perhaps with different assumptions about the fre- 
quency dependence of the amplitudes; they are urgently 
needed to strengthen the case for the reality of this signal. 
However, it is certainly encouraging that evidence for g 
modes appears to be present in the GOLF data, and fur- 



ther analyses and simulations in this direction will be of 
great interest. 



5. NEAR-SURFACE EFFECTS 



The near-surface regions of the Sun are obviously highly 
complex, with a strong interplay between thermal struc- 
ture, flows and magnetic fields. The richness and chal- 
lenges of this region are reviewed by Erdelyi (these pro- 
ceedings). Thus the ability to make detailed investiga- 
tions of these regions by local helioseismology, reviewed 
by Toomre (these proceedings) and discussed by a num- 
ber of papers in the present volume, is extremely inter- 
esting. At the photospheric level, and presumably be- 
low, magnetic effects are strongly localized to active re- 
gions where the magnetic and gas pressure are compara- 
ble. Particularly important is the critical surface where 
the sound speed equals the Alfven speed and where con- 
sequently mode conversion occurs. The detailed interac- 
tion between the magnetic field and the acoustic waves 
depends critically on the inclination angle of the field 
and the direction of propagation of the waves 19211 . Jef- 
feries et al. (these proceedings) have found that such ef- 
fects allow the propagation of waves into the solar chro- 
mosphere at frequencies below the acoustical cut-off fre- 
quency (around 5 mHz) in the non-magnetic part of the 
atmosphere. Such 'magneto-acoustical portals' are asso- 
ciated with the small-scale strong magnetic features that 
are ubiquitous on the solar surface. Jefferies et al. demon- 
strate that the wave energy transmitted through these re- 
gions at relatively low frequency may play a substantial 
role in the heating of the chromosphere. Furthermore, 
this effect likely has a significant influence on the damp- 
ing, and possibly the frequencies, of global modes of os- 
cillation. 

A striking parallel to these solar effects, albeit on a far 
larger scale, is discussed by Kurtz et al. (these proceed- 
ings), in connection with the rapidly oscillating Ap stars. 
Here the properties of the global modes are dominated 
by the large-scale magnetic field in these stars. Strik- 
ingly, the strong separation of elements with altitude, and 
across the stellar surface, allows detailed inferences to 
be made of the variation of the oscillation eigenfunctions 
with height in the atmosphere, thus providing potentially 
very stringent tests of the understanding of the physics of 
the modes under the influence of the magnetic field. 

The interaction between the magnetic fields and the 
acoustic waves poses particular problems for the inter- 
pretation of results of local helioseismology in the vicin- 
ity of active regions. Lindsey & Braun 1 70] identified an 
acoustical 'shower-glass effect', resulting from the dis- 
tortion of the acoustic waves by subphotospheric mag- 
netic fields, which impairs the analysis, at least in the case 
of helioseismic holography (see also Lindsey, these pro- 
ceedings). These effects are further discussed by Schun- 
ker & Cally (these proceedings), while Zhao (these pro- 
ceedings) notes that similar effects are present in time- 
distance analyses, although possibly less severe. Related 



problems concern the interpretation of results of local 
helioseismology in terms of direct magnetic, or induced 
thermal, effects. These issues must be understood before 
a fully reliable interpretation, in physical terms, can be 
made of the results of local helioseismology. 

The goal of helioseismology, global as well as local, is 
obviously to resolve the properties of the solar interior. 
However, the physical properties of the observed modes 
unavoidably place constraints on the resolution that can 
be achieved. Unresolved aspects of the near-surface ef- 
fects must somehow be suppressed in the analysis of the 
observations to infer properties deeper in the solar inte- 
rior. In the case of structure inversions in global helio- 
seismology it is well known that the dominant difference 
between the observed and the model frequencies is gen- 
erally caused by unresolved effects near the surface (e.g., 
ll3~fll ). From the frequency dependence of these effects 
they must predominantly be localized to the uppermost 
parts of the solar interior and the solar atmosphere, with 
likely culprits being convective effects where the strati- 
fication is strongly superadiabatic or effects of the mode 
damping and excitation. In particular, modes of relatively 
low frequency, which are reflected below this region, suf- 
fer minimal effects. These effects are unlikely to be fully 
resolved in the radial direction and hence must be sup- 
pressed in the inverse analysis; in the case of global helio- 
seismology this is largely possible owing to the availabil- 
ity of modes over a broad range of frequency and degree. 
The frequency effects can be regarded as corresponding 
to a phase shift of the waves when they are reflected near 
the surface; as such they should be present also in the 
properties of waves as studied in local helioseismology, 
including time-distance analysis; it is interesting that this 
has recently been confirmed 1 2311 . 

These near-surface effects pose serious problems in the 
analysis of data on solar-like oscillations in distant stars, 
where only data on low-degree modes are available. Even 
here the effects, insofar as the investigations of stellar 
cores are concerned, can to some extent be su ppre ssed 
through a suitable combination of the data I85l 1771 18411 . 
However, a great deal of information, particularly about 
the global properties of the star such as its radius, is 
lost in the process. Based on the solar experience, one 
may expect that lower-frequency modes are less affected 
and hence provide largely unbiased diagnostics; however, 
these modes generally have lower amplitudes and hence 
are challenging to observe. 

Evidently, the ultimate goal is to understand these effects 
and hence use them as additional diagnostics. Promis- 
ing progress in this regard has been made on the basis 
of hydrody namical simulations of convection in the so- 
lar 1 82, 81] and stellar Ill02ll case. The results presented 
by Stein et al. (these proceedings) have taken such sim- 
ulations to a new level of detail. On this basis it should 
be possible to investigate directly the effects of convec- 
tion on the properties of the modes, by analysing simu- 
lations extending over a sufficiently long period that the 
modes excited in the simulation has time to reach a suit- 
able equilibrium. This could obviously be done for a 
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Figure 3. Relative differences in squared sound speed, in 
the sense ( Sun ) - ( model), against fractional distance r/R 
to the solar centre. The open diamonds are for Model S 
h320 . The filled diamonds are for a model with the revised 
composition presented by Grevesse, while the crosses are 
for a model where in addition the neon abundance has 
been increased by 0.5 dex. The (barely visible) verti- 
cal bars show 1 — a errors, whereas the horizontal bars 
give a measure of resolution in the underlying inversion. 
(Houdek et al., in preparation.) 

range of stellar parameters, thus possibly allowing a reli- 
able calibration of simpler formulations for more general 
use. Also, additional effects such as magnetic fields could 
be included in the simulations. 



6. COMPOSITIONAL CHALLENGES 



Helioseismology has provided what must be regarded as a 
relatively accurate determination of properties of the so- 
lar internal structure, particularly the sound speed. Al- 
though the inferences most often are made with reference 
to solar models, leading to a determination of the dif- 
ferences in, e.g., sound speed between the Sun and the 
model, the final sound-speed profile resulting from ap- 
plying that correction to the model is largely insensitive 
to the reference model fun . There remain issues of res- 
olution and error correlation which mean that testing of 
a specific model is best carried out by using that model 
as reference in an inverse analysis to determine the dif- 
ferences with the solar structure, as discussed in m. 
However, in general helioseismology has been remark- 
ably successful in determining the solar internal sound 
speed. 

The success of a solar model then depends on how well 
it reproduces that helioseismic determination. In paral- 
lel with the development of helioseismology calculations 
of solar evolution were improved, leading to models that 
were in reasonable agreement with helioseismic infer- 
ences, although still differing from them by much more 
than the formal errors in the inferences (see Figure|3j; a 
typical example is 'Model S' I32T which has been exten- 



sively used in helioseismic analyses. Although the im- 
provements in the computations were inspired by helio- 
seismology no explicit parameter adjustments were made 
to match the helioseismic results; the input to the calcu- 
lation was obtained from the, at the time, best possible 
knowledge about the physics of the solar interior and the 
other observables, such as the ratio Z s /X s between the 
present surface heavy-element and hydrogen abundances. 
By adjusting parameters or aspects of the internal physics 
models in even closer agreement wit h the helioseismic 
sound speed can be constructed (e.g., Il06l0 . 

Recently, a major revision of the determinations of solar 
surface abundances has taken place 0, IH l53ll . as sum- 
marized by Nicolas Grevesse at the conference. The 
improvements in the abundance analyses involve using 
three-dimensional and time-dependent hydrodynamical 
simulations with a realistic treatment of radiative trans- 
fer, as well as abandoning the assumption of local ther- 
modynamic equilibrium in the treatment of the formation 
of some spectral lines. A strong indication of the im- 
provement in the line-formation modelling resulting from 
the use of the hydrodynamical models is evident from 
the fact that the line profiles can be modelled without 
the use of ad hoc micro and macro turbulence. Also, 
strong support for these new results come from the fact 
that, unlike earlier determinations, consistent abundances 
are obtained for a given element from different spectral 
lines. The most dramatic change is in the abundance of 
oxygen which is reduced by a factor of 0.54 in the new 
determination, relative to the values used until recently 
in solar modelling |2]; also, the carbon abundance is re- 
duced by a factor of 0.68. As a result Z s /X s is changed 
from the value 0.0245 used, for example, in Model S to 
Z s /X s = 0.0165. 

The heavy-element abundances predominantly affect so- 
lar structure through the opacity; in fact, oxygen makes 
a dominant contribution to the opacity just beneath the 
convection zone. Thus it is not surprising that the change 
in composition leads to a change in solar structure, as 
illustrated in terms of sound speed in Figure |3j clearly, 
the composition change has caused a strong deterioration 
in the agreement between the model and the Sun. Also, 
the initial hydrogen abundance required to calibrate the 
model is changed, leading to a decrease in the helium 
abundance Y c in the convective envelope from 0.245 in 
Model S to 0.221, and hence to a present surface heavy- 
element abundance Z s = 0.0125 (compared with 0.0181 
in Model S). Finally, the depth d cz of the convective en- 
velope changes from O.289i?0 (Rq being the solar ra- 
dius) in Model S to O.271i?0. As discussed, for exam- 
ple, in lll4Tl the values for Y e and d cz obtained for models 
like Model S were consistent with helioseismic inference 
while the new values are not. 

This renewed discrepancy between solar models and the 
helioseismically determined solar structure is a serious 
problem for solar and stellar modelling. Guzik (these pro- 
ceedings) discusses it in more detail, including the many 
attempts that have been made to solve it. There is clearly 
a need for independent confirmation of the analysis by 



Asplund et al.; in fact, the result has been questioned 
on the basis of possible problems with the solar atmo- 
spheric temperature profile used 1 10]. On the assumption 
that the revised abundances are correct, it is evident that, 
since the problem arises from the change in opacities re- 
sulting from the change in composition, it can be solved 
by an intrinsic opacity modification that compensates for 
the change in the heavy-element abundance. The mag- 
nitude of the required opacity change is probably larger 
than the quoted precision of the opacity calculations, al- 
though perhaps not larger than the real uncertainty in the 
computed opacities. Such a change is hardly satisfying 
without a physical justification; however Kurucz (private 
communication) has noted that the very many weak lines 
that are neglected in current opacity calculations might 
add up to a contribution of the required size, thus strongly 
motivating consideration of such effects. Opacity may 
also be increased by increasing abundances of other ele- 
ments making substantial contributions to the opacity. In- 
teresting examples are neon and argon, whose solar abun- 
dances are very uncertain, since they have no spectral 
lines formed in the solar photosphere. It has been shown 
that, for example, an increase in the neon abundance by 
0.5 dex (i.e., a factor 3.16) comes close to restoring the 
agreement between the model and the Sun seen for Model 
S ifl- lol : an example of such a model is also illustrated in 
Figure[3] However, although there may be evidence for a 
such a higher abundance of neon in solar-neighbourhood 
stars 1 39] the applicability of these results to the Sun has 
been questioned (e.g., I89B '). A further test of the abun- 
dance of neon in the solar convection zone can in prin- 
ciple be obtained from helioseismic constraints on the 
effect on the equation of state, particularly the adiabatic 
compressibility 1 69] . An initial analysis indicates that the 
required enhanced neon abundance is inconsistent with 
the helioseismic results @|. 

Guzik (these proceedings) considers several other poten- 
tial solutions to the problem posed by the new abun- 
dances, including energy transport against the tempera- 
ture gradient by gravity waves, hence corresponding to a 
higher effective opacity (interestingly, gravity waves have 
also been invoked to account for the transport of angu- 
lar momentum; see Talon, these proceedings). She con- 
cludes that the previous reasonable (though not perfect) 
agreement between the solar models and the Sun might 
have arisen from a fortuitous cancellation of errors, and 
that the new discrepancies are a good motivation to re- 
consider the computation of solar models, including also 
normally neglected hydrodynamical effects. This clearly 
represents an exciting, if difficult, challenge. 



7. PHYSICS OF STELLAR INTERIORS 



Opacity is not the only aspect of the 'microphysics' used 
in stellar modelling that might be in question. The is- 
sue of electron screening of nuclear reactions in stel- 
lar interiors has seen vi goro u s dis cussion over an ex- 
tended period (e.g., 0, lH, OH, El)- Most calcu- 
lations of stellar evolution use the classical description of 



the screening enhancement developed by Salpeter |86], 
based on an equilibrium treatment of the average charge 
around the reacting nuclei. This description may, how- 
ever, reasonably be questioned on the basis of the relevant 
timescales and the moderate number of particles typically 
involved in the screening. This has motivated investi- 
gations of the interactions through molecular-dynamics 
calculations, particularly by Shaviv and Shaviv 1 961 19511 . 
Given the complexity of such calculations, it is very en- 
couraging that Mussack et al. (these proceedings) have 
embarked on similar but entirely independent calcula- 
tions. The results are eagerly awaited. 

Although the microphysics is complex and challeng- 
ing in stellar evolution calculations, even more serious 
challenges are caused by potential hydrodynamical and 
magneto-hydrodynamical phenomena in stellar interiors. 
Convection is an obvious difficulty; as discussed by Stein 
et al. (these proceedings) an apparently realistic descrip- 
tion is possible of the near-surface convection, although 
at large computational expense. Also, Toomre (these pro- 
ceedings) presents results of large-scale simulations of 
convection and rotation in the solar convection zone, in- 
cluding also magnetic effects, although not yet giving a 
fully adequate modelling of, for example, the solar inter- 
nal differential rotation. Such detailed simulations of the 
deeper interiors of stars are hampered by the vast range of 
timescales that need to be considered, between the over- 
all thermal timescale and the dynamical timescales. Thus 
no fully realistic simulations have been made, for exam- 
ple, of convective penetration beneath the solar convec- 
tion zone. 

Observationally, it is clear that non-standard phenomena 
must be considered in stellar modelling. The nearly uni- 
form rotation of the solar radiative interior requires mech- 
anisms for the transport of angular momentum towards 
the surface, to reduce internal rotation from the assumed 
initial rapid rate in the early Sun. Also, in many cases the 
observed surface abundances of stars require mixing be- 
yond the convectively unstable regions, hence indicating 
some kind of hydrodynamical instability. Palacios et al. 
(these proceedings) provide an overview of the required 
non-standard extensions of stellar modelling. 

A dominant source of dynamical effects in many stars 
is rotation. It is relatively straightforward to include the 
spherically symmetric component of the centrifugal force 
in stellar modelling, although this requires a model for 
the evolution of the internal rotation rate with age. A 
much more difficult issue is the meridional circulation 
and the instabilities which are likely caused by rotation 
and which lead to transport of chemical elements and an- 
gular momentum. Zahn 1 1 10] developed an approximate 
treatment of these processes which has been refined and 
extensively applied. Interestingly, the resulting transport 
of angular momentum is insufficient to account for the 
solar internal rotation (see Palacios, these proceedings). 

Schatzman |87] proposed that gravity waves excited by 
penetration from the solar convective envelope might be 
efficient in transporting angular momentum in the radia- 



tive interior, thus accounting for the internal solar rota- 
tion. As discussed by Talon (these proceedings) further 
refinements, including differential filtering of the waves 
in a layer just below the convection zone, are required 
to make this a potentially viable mechanism. The inter- 
nal rotation rate resulting from such calculations is quite 
similar to the helioseismically inferred behaviour. Alter- 
natively, magnetic fields may play an important role in 
coupling the interior to the convection zone. Gough & 
Mclntyre 15211 pointed out that even a weak primordial 
field would be sufficient and furthermore discussed how 
this could account for the localized transition of rotation 
in the solar tachocline (see Toomre, these proceedings). 
Alternatively, the field may be generated by dynamo ac- 
tion through differential rotation in the radiative interior 
119911 ; such a model also results in near uniform rotation 
of the radiative interior, as observed in the Sun H43I1 - 

The processes discussed here should of course be in- 
cluded in a complete modelling of stellar structure and 
evolution. Mathis et al. (these proceedings) describe 
an ambitious programme to develop a stellar evolution 
code making this possible, with an initial application to 
a 1.5 Mq star. Development and verification of such a 
code is evidently a major task which is urgently needed, 
given the extensive asteroseismic results expected from 
upcoming missions. The observational tests of the re- 
sults will require analysis of a broad range of observables 
in different types of stars. The solar internal structure 
and rotation obviously provide strong constraints. Addi- 
tional constraints are available from the detailed distribu- 
tion of element and isotope abundances observed in dif- 
ferent stars, reflecting the internal nuclear reactions and 
the mixing processes bringing the products of these reac- 
tions to the surface, as well as the destruction of lithium 
and possibly beryllium through nuclear burning. Also, 
the detailed asteroseismic results expected from CoRoT 
and Kepler, and eventually from network projects such 
as SONG (see Section II lb . should provide strong con- 
straints on the modelling. 

It should be kept in mind that the current and planned 
level of modelling, although relatively complex, still re- 
lies on highly simplified descriptions of the magneto- 
hydrodynamical processes in stars. Thus the modelling 
should be complemented by investigations aiming at im- 
proving our understanding of these processes, to improve 
their description in the stellar evolution codes. This will 
undoubtedly involve substantial numerical simulations, 
the results of which, while perhaps not directly applicable 
to stellar conditions, can at least be extrapolated to such 
conditions. 



8. MODE EXCITATION AND DAMPING 



Two different mechanisms are available to excite a given 
mode: self-excitation where the mode is linearly unsta- 
ble, with an amplitude that consequently grows as a re- 
sult of the conversion of thermal energy into mechanical 
energy; and external forcing of an intrinsically damped 



mode. In the former case a critical region of the star acts 
as a heat engine, with the appropriate phasing of heat- 
ing and compression during the pulsation cycle. In or- 
der for the effect to lead to overall instability that region 
must be located at a suitable depth within the star; con- 
sequently, stars pulsating because of heat-engine excita- 
tion tend to be located in relatively well-defined regions 
of the HR diagram where that condition is satisfied. A 
typical example is the classical Cepheid instability strip 
where the excitation is caused by opacity features associ- 
ated with the second ionization of helium. The external 
forcing of stable modes can take place through tidal in- 
teraction in binary systems. In single stars, which are 
the focus here, the only mechanism that has been consid- 
ered in any detail is forcing by near-surface convection; 
it should be recalled, however, that excitation of gravity 
waves through convective penetration beneath convective 
envelopes may play a role for stellar structure and evolu- 
tion (cf. Section[7}. 

Dupret et al. (these proceedings) review the properties of 
modes excited through the heat-engine mechanism. In 
relatively hot stars, where convection at most makes a 
minor contribution to energy transport in the relevant re- 
gions of the star, the critical aspect of the energy transport 
leading to excitation is the behaviour of opacity, and one 
therefore often talks about the 'k mechanism' . Excitation 
is typically associated with rapid variations, or 'bumps', 
in the opacity near a given temperature; instability arises 
for modes such that the oscillation period matches the 
thermal timescale of that part of the star which is out- 
side the region where local excitation takes place OH- 
Each characteristic bump in the opacity is generally asso- 
ciated with one or more classes of pulsating stars; longer- 
period modes are evidently preferred when the bump oc- 
curs deeper in the star, i.e., in cooler stars. As discussed 
below, examples of this behaviour are found in the region 
of the (3 Cephei and slowly pulsating B stars, as well as 
amongst the pulsating subdwarf B stars. 

In modelling cooler pulsators, the treatment of the inter- 
action between convection and pulsations is a major chal- 
lenge. This concerns both the perturbation to the con- 
vective flux, which affects the heating, and the pertur- 
bation to the turbulent pressure, which has a direct dy- 
namical effect on the pulsations. As discussed by Dupret 
et al. convective effects are likely responsible for the re- 
turn to stability at the cool edge of the Cepheid instability 
strip. This is confirmed by a number of independent cal- 
culations, using rather different physical descriptions of 
how pulsations affect convection. However, the detailed 
cause of the damping differs between the different formu- 
lations; in some calculations (e.g., |58]) the perturbation 
to the turbulent pressure dominates the damping, while 
in others (e.g., [41]) damping is dominated by the con- 
vective flux. It is evidently essential that these issues be 
clarified through detailed comparisons between the dif- 
ferent calculations. Also, it is important to ensure that 
the convection is treated in a manner that is consistent 
between the calculation of the equilibrium model and the 
pulsations. 



In the 7 Dor stars, beyond the cool side of the instabil- 
ity strip, the mode excitation is apparently caused by a 
heat-engine mechanism, but controlled by the effect of 
convection. In the critical layer at the base of the con- 
vective envelope the convective timescale is much longer 
than the pulsation period, and hence convection does not 
react to the pulsations. This 'frozen convective flux' leads 
to heating just beneath the convection zone in the appro- 
priate phase of the pulsations, and hence to instability 
I56I1 . Dupret et al. note that calculations with a detailed 
treatment of the reaction of convection confirm this sim- 
ple picture. It should be noticed that this is an example 
of heat-engine excitation that cannot be regarded as an 
opacity mechanism. 

Heat-engine excitation demonstrates instability of a mode 
but says nothing about the limiting amplitude. Except for 
large-amplitude pulsators, where the excitation mecha- 
nism apparently saturates, there is no definite understand- 
ing of the amplitude-limiting mechanism or the resulting 
distribution of mode amplitudes. In a detailed analysis 
Nowakowski 17511 showed that resonant mode interaction 
between unstable and stable modes was not able to ac- 
count for the observed amplitudes, as had previously been 
suggested. Thus a reliable theory for the mode selection 
in heat-engine pulsators remains an important challenge 
to the study of such stars. 

The stochastic excitation of solar-like oscillations is re- 
viewed by Houdek (these proceedings). In this case the 
amplitudes are determined by the balance between the en- 
ergy input from convection and the damping rate, leading 
to a definite prediction of the pulsation amplitudes. The 
calculation of the damping rate evidently requires treat- 
ment of convective effects which, as already mentioned, 
is uncertain. In the solar case, however, accurate obser- 
vational determinations of the damping rates have been 
obtained from the measured line widths of the modes. As 
illustrated by Houdek, computed damping rates, based on 
the convection theory of Gough as further developed by 
Balmforth 15(1 11511 - can be adjusted to be in reasonable 
agreement with these observations. This gives some hope 
that the description may be extended to other stars. In 
fact, the computed damping rates for a Cen A and B are 
similar to those inferred from observations (and to the so- 
lar values). On the other hand, for the giant star £ Hydrae 
the computed lifetime 15911 is considerably longer than 
the lifetime i nferr ed, albeit somewhat indirectly, from 
observations 1 1 111 . The observational determination of 
lifetimes will be greatly improved by the expected long- 
duration observations from, e.g., the CoRoT and Kepler 
missions (see Section fTTl . which will allow a direct de- 
termination of the line widths. However, it is obvious that 
improvements are also required in the computation of the 
damping rates. 

The energy input comes partly from dynamical effects 
through the Reynolds stresses, and partly from thermo- 
dynamical effects through fluctuations in the buoyancy 
force (often referred to as entropy fluctuations). Compu- 
tation of these contributions obviously requires a model 
for the temporal and spatial spectra of convection in 
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Figure 4. Energy input in stochastic excitation of solar 
acoustic modes as a function of depth below the photo- 
sphere and cyclic frequency v. The upper panel shows 
results for a simple mixing-length model of convection, 
on a linear scale, whereas the lower panel, on a loga- 
rithmic scale, is based on hydro dynamical simulations of 
near-surface convection [100]. (Adapted from Houdek, 
these proceedings.) 

the equilibrium model, as discussed by Houdek. Since 
this involves higher-order moments of the turbulent field, 
appropriate closure models are also required; an inter- 
esting recent development is the use of closure models 
that explicitly take into account the asymmetry between 
broad upwellings and narrow downdrafts in stellar con- 
vection [21]. Further insight can be obtained from com- 
paring the simple models, typically based on some ver- 
sion of mixing-length theory, with the energy input com- 
puted from detailed hydrodynamical models of convec- 
tion UPO], as illustrated in Figure^ It is evident that the 
general behaviour is similar in the two cases, although ex- 
citation in the simulations extends to considerably greater 
depth than in the simple model. 

In the solar case, the computed amplitudes are generally 
in reasonable agreement with observations; however, if 
the computed damping rates are used the predicted am- 
plitudes are substantially larger than the observed values 
at low frequencies, where in fact the computed damping 
rates are too small. Using instead the observed damp- 
ing rates, the predicted and observed amplitudes are in 
good agreement |28[. For other stars the situation is less 
clear, as summarized by Houdek. In particular, it appears 



that the predicted amplitudes are substantially higher than 
observed for stars hotter than the Sun, such as Procyon. 
Interestingly, Houdek notes from hydrodynamical simu- 
lations that there may be some cancellation, so far not 
included in the calculations, between the contributions 
from Reynolds stress and buoyancy to the excitation; this 
can perhaps explain at least partly why the computed val- 
ues are too high. In the case of £ Hydrae it is striking that 
the predicted amplitude is in excellent agreement with 
observations, despite the large discrepancy in the damp- 
ing rates. 

An issue closely related to mode excitation, also consid- 
ered by Dupret et al. and by Houdek, is the relation be- 
tween observable properties of the oscillations as seen in 
intensity or radial velocity. This is important for the use 
of different observations of a given mode to identify the 
degree of the mode. The observed properties obviously 
depend on the response of the atmosphere to the oscil- 
lations. At frequencies well below the acoustical cut-off 
frequency a quasi-static treatment of the oscillations in 
the atmosphere is probably adequate 1 40], while at higher 
frequencies, such as are observed in the Sun and solar- 
like pulsators, a fully time-dependent calculation may be 
required Il29ll . Also, the treatment of the convective flux 
is obviously important in cooler stars; an interesting anal- 
ysis of the effect on mode identification in the S Scuti star 
FG Vir was carried out in [36]. Additional complications 
in the mode identification arise in rapidly rotating stars, 
where a given mode may have contributions from more 
than one spherical-harmonic degree 13511 . 

It is evident that the near-surface and atmospheric re- 
gions are extremely complex in stars with convective en- 
velopes, yet crucial for the understanding and observa- 
tions of oscillations in such stars; this also, as discussed 
in Section |5] applies to the interpretation of the oscilla- 
tion frequencies. Although the theoretical understanding 
of these regions is improving, there is no doubt that ma- 
jor advances will result from the hydrodynamical simu- 
lations. The new simulations of solar near-surface con- 
vection reported by Stein et al. (these proceedings) are a 
gold-mine for the further development of studies of ex- 
citation and other properties of stellar oscillations. Os- 
cillations are spontaneously excited in the simulation and 
the duration, 65 hours of solar time, may be sufficient to 
allow a determination of the damping rate of the modes 
and an understanding of the dominant factors controlling 
the excitation and damping II49I1 . With similar simula- 
tions for selected examples of other stars we can hope to 
be able to calibrate simpler models and in this way obtain 
reliable calculations of mode properties for a wide range 
of stars. When compared with the upcoming large ob- 
serving campaigns this provides excellent prospects for 
an improved understanding of the properties of the oscil- 
lations and, one may hope, observational constraints on 
the properties of convection. 



9. PROGRESS ON ASTEROSEISMOLOGY 



The remarkable recent development in asteroseismology 
covers most of the HR diagram and a broad range of stel- 
lar evolutionary stages. Developments in ground-based 
techniques for radial-velocity observations have con- 
tributed greatly to this, as have substantial programmes 
for coordinated observation of selected objects. Also, 
as reported by Jaym ie Matthews at the conference the 
MOST satellite 1 109] has produced striking results for a 
substantial number of stars, although perhaps not quite 
reaching the expected level of sensitivity towards solar- 
like oscillations. 

Bedding & Kjeldsen (these proceedings) discuss the 
stochastically excited solar-like oscillations which are 
likely found in all stars with extensive outer convection 
zones. It is striking that such oscillations, with the char- 
acteristic Lorentzian profiles, have likely been detected 
in red supergiants such as Betelgeuze (a Ori), from ama- 
teur observations recorded by the American Association 
of Variable Star Observers (see also |65]). Solar-like os- 
cillations in red giants are well established; a remarkable 
set of data from the MOST satellite has been obtained 
for the star e Ophi (Barban et al., these proceedings). In- 
terestingly, Hekker et al. (these proceedings) find from 
radial-velocity observations of the same star that the ob- 
served modes are likely non-radial; this raises problems 
for the theoretical understanding of oscillations in red gi- 
ants, according to which such modes are damped substan- 
tially more strongly than are radial modes 14211 . 

Closer to the main sequence detailed frequency spectra 
have been obtained for a substantial number of stars. An 
interesting case is /i Ara which also hosts a planetary 
system and, like many other planetary hosts, has a rel- 
atively high heavy-element abundance Z. As discussed 
by Vauclair (these proceedings) the asteroseismic anal- 
ysis in principle allows to distinguish between models 
with a globally high Z and models where just the con- 
vection zone is so enriched; in practice the current data, 
although extensive, are still not quite sufficiently detailed 
to allow a decision on this issue (see also 1 20] ). In an 
analysis of an extended series of data from the WIRE 
satellite Fletcher et al. (these proceedings; see also 1 44]) 
were able to determine directly the average mode lifetime 
and rotational splitting for a Cen A. Potentially very in- 
teresting, although controversial, results on ?; Boo have 
been reported based on MOST observations 1 55]; these 
indicate a sequence of radial modes extending from the 
previously observed modes towards modes of low or- 
der. Such modes would have very valuable diagnostic 
potential; however, their detection at relatively high am- 
plitudes is contrary to our understanding of stochastic ex- 
citation and entirely unlike the properties of the spectra 
observed in other stars with solar-like oscillations. Thus 
some scepticism towards this claimed detection is proba- 
bly warranted. 



In contrast to the solar-like oscillators, the heat-engine 
pulsators (see also Section |8j occupy relatively well- 



defined instability regions in the HR diagram, determined 
by the physical mechanisms responsible for the excita- 
tion. Also, these stars show pulsations at larger ampli- 
tudes than the solar-like pulsators near the main sequence 
(however, smaller-amplitude pulsations, below the cur- 
rent observational threshold, are likely also present; e.g., 
112410 . Matthews presented extensive results from MOST 
for a broad range of stars, including 6 Scuti and (3 Cephei 
stars; very interesting results have been obtained for the 
rapidly oscillating Ap star HD 1217, showing additional 
'anomalous' frequencies affected by magnetic effects. 
Aerts (these proceedings) discusses oscillations in hot 
main-sequence stars. Detailed analyses of (3 Cephei stars 
have become possible thanks to very extended sets of co- 
ordinated observations over the last few years. An inter- 
esting case is v Eri where standard models do not find 
instability of all the observed modes in models that fit 
the observed frequencies l9 H78ll : the instability is likely 
caused by enhancements of the abundances of iron-group 
elements in the driving region, through settling and ra- 
diative levitation, as is also invoked for the subdwarf B 
pulsators (see below). Significant asteroseismic results 
can be obtained from even a small number of modes, if 
these are properly identified; thus Mazumdar et al. 17211 
obtained stringent constraints, including a determination 
of the core overshoot, for the (3 Cephei star (3 CMa from 
just three modes of which two were shown to be undergo- 
ing an avoided crossing. Slowly pulsating B (SPB) stars 
have been discovered in substantial numbers but their 
very long periods (of the order of a day or more) make 
a detailed analysis very difficult and time consuming. 
Also, while the identification of the observed frequencies 
with high-order g modes makes them potentially interest- 
ing diagnostics, the very dense spectrum of theoretical 
frequencies and the generally rapid rotation (compara- 
ble to the oscillation frequencies) of the stars complicate 
the interpretation of the observations. An extensive set 
of frequencies in the star HD 163830 has been obtained 
from the MOST observations, but a definite identification 
of the degree of the modes, let alone identification with 
modes in stellar models, has not yet been possible 01 • 

Fontaine (these proceedings) discusses the pulsating sub- 
dwarf B (sdB) stars, a recently identified class of variable 
stars. These are helium-burning stars at the extreme blue 
end of the horizontal branch who have lost almost all of 
the hydrogen-rich envelope, possibly as a result of binary 
evolution. The generally extensive sets of observed fre- 
quencies provide an excellent potential for studying the 
properties of these highly evolved stars and hence con- 
straining the earlier evolutionary phases. Two classes of 
pulsating sdB stars have been found: short-period vari- 
ables pulsating in p modes and, at slightly lower effec- 
tive temperature, long-period variables pulsating in high- 
order g modes. The excitation of the modes can be under- 
stood in terms of heat-engine driving in the opacity bump 
associated with iron-group elements, but only if the abun- 
dances of these elements are enhanced by radiative levi- 
tation 1 45]. The required enhancement is established on a 
timescale short compared with the evolutionary timescale 
of these stars 



10. ASTEROSEISMIC DIAGNOSTICS 



Reliable extraction of information about stellar proper- 
ties from the observed oscillation frequencies requires a 
good understanding of the relation between those prop- 
erties and the properties of the oscillation frequencies. 
This greatly aids the combination of frequencies in ways 
that isolate measures of specific aspects of the stellar in- 
teriors, with an understanding of how the observational 
uncertainties give rise to uncertainties in these measures. 
An obvious example is the extensive use of the large and 
small frequency separations in asteroseismic diagnostics, 
based on solar-like oscillations, to determine global pa- 
rameters of the stars. 

Evidently the goals of asteroseismology, particularly 
given the upcoming observational projects, go much be- 
yond such basic global parameters. Consequently, we 
need to understand how the finer details of the observed 
frequencies are related to more subtle aspects of the stel- 
lar interior. An important example, discussed by Houdek 
& Gough (these proceedings), is the effects on the fre- 
quencies of relatively sharp features, denoted 'acoustical 
glitches' by Houdek & Gough. Such features include the 
effect on the adiabatic compressibility of the ionization 
of helium, and the abrupt change in the temperature and 
hence sound-speed gradient at the base of a convective 
envelope. Houdek & Gough develop analytical expres- 
sions for the influence of these effects on the second dif- 
ference of frequencies of low-degree modes of like de- 
gree. The expressions contain parameters that, for exam- 
ple, provide information about the helium content of the 
convective envelope and its depth. They provide excel- 
lent fits to second differences of computed frequencies 
for a solar model. The quality of the fit, and the phys- 
ical basis for the parametrization, is undoubtedly supe- 
rior to previous examples of similar analyses Ii73lfl9i 17411 . 
Thus it seems likely that the parametrization by Houdek 
& Gough will provide a more reliable determination of 
these properties in stellar interiors, once adequate data 
are available. 

Rigorous results in the theory of stellar oscillations are 
rare, for general stellar models. In most cases insights 
into the properties of the modes are based on asymptotic 
analysis or numerical calculations. Thus the announce- 
ment by Takata 1 103] of an identity strictly satisfied be- 
tween the perturbation quantities of oscillations of degree 
/ = 1 is a major advance in the theory of stellar oscil- 
lations. Such dipolar oscillations are special in having 
a geometry where the centre of mass of each perturbed 
shell in the star is displaced, while obviously the centre 
of mass of the star as a whole remains stationary dur- 
ing the oscillation. This property forms the basis for the 
identity. As a result of the identity, the equations of adi- 
abatic oscillations can be reduced from the usual fourth- 
order system to a second-order system. In these proceed- 
ings Takata further develop this theory in a conceptually 
simpler form, which can be directly compared with the 
classical asymptotic description in the Cowling approxi- 
mation of neglecting the perturbation to the gravitational 
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Figure 5. Computed dimensionless frequencies a for 
dipolar modes for a IMq evolution sequence against 
central condensation; here a = (i? 3 /GM) 1 / 2 cj, d be- 
ing the angular frequency, R and M the stellar radius 
and mass and G the gravitational constant, and central 
condensation is measured as the ratio between the central 
and average densities p c and p av . Small dots mark modes 
that are labelled correctly by the usual scheme i9Jl l7dl/ . 
whereas crosses mark modes where this scheme fails. The 
dashed line is defined by a 2 — p c /p av . (Christensen- 
Dalsgaard & Takata, in preparation.) 

potential. This will undoubtedly improve our understand- 
ing of the diagnostic potential of these modes; it is worth 
recalling that amongst acoustic modes the dipolar modes, 
together with the radial modes, are those that penetrate 
most deeply into the stellar core. 

The reduction of the order of the equations to two also 
allows a simple procedure for defining the radial order of 
the dipolar modes. The commonly used procedure, es- 
sentially based on assuming the Cowling approximation, 
fails for models of high central condensation (e.g., l54ll . 
IBOIO . In contrast, applying a similar procedure to the 
variables satisfying the exact second-order equations de- 
rived by Takata allows a robust assignment of order for 
any model. Such a secure labelling further helps the un- 
derstanding of how the modes change as the star evolves 
and allows tests to ensure that no modes have been missed 
in computations of stellar oscillation spectra. As an ex- 
ample, Figure[5]illustrates where the previously used pro- 
cedure fails in an evolution sequence for a 1M@ model. 



11. CHALLENGES AND PROSPECTS 



It is obvious that great progress has been made to- 
wards an improved understanding of the interior of the 
Sun and other stars, through extensive observations and 
modelling. Equally obviously, many challenges remain. 



The treatment of rotation remains a critical issue, de- 
spite the progress on including rotational effects in stel- 
lar modelling. For the many stars rotating substantially 
faster than the Sun the treatment of the effect of rota- 
tion on the oscillation frequencies requires going beyond 
the first-order expansion typically used in helioseismol- 
ogy. An expansion up to third order has been devel- 
oped 1 98] but has so far, perhaps because of its complex- 
ity, seen relatively limited use. Also, even this expan- 
sion is inadequate for many of the stars for which aster- 
oseismic investigations are carried out; here a fully two- 
dimensional calculation is required to match the obser- 
vational precision 1 68, 80]. Similarly, a consistent fully 
two-dimensional calculation of stellar structure and evo- 
lution will be needed; some prog ress has been made in 
this direction (e.g., 1 38], |83], |63]), although much work 
obviously remains before this can be combined with the 
treatment of hydrodynamical effects discussed in Sec- 
tion The development of this required next level of 
stellar modelling will undoubtedly take place through a 
combination of idealized, if highly complex, hydrody- 
namical simulations addressing specific aspects of the 
physics and comprehensive attempts to model the full 
stellar evolution with somewhat simplified representation 
of the physics. 

The observational prospects for helio- and asteroseismol- 
ogy are bright. As argued in Section |2] the BiSON and 
GONG networks should be kept operating to ensure co- 
herent data covering several solar cycles. With the ex- 
tension of SOHO operations through 2009 a reasonable 
overlap should be ensured with the observations from the 
Solar Dynamics Observatory (SDO). The launch of SDO, 
expected in late 2008, will result in a dramatic increase in 
data on solar waves and oscillations, with Helioseismic 
and Magnetic Imager (HMI) yielding full-disk Doppler- 
velocity data with a spatial resolution of 1.5 arc seconds 
and a time resolution of less than one minute. Unlike the 
MDI instrument on SOHO such high-resolution data will 
be available continuously throughout the mission. This 
will revolutionize the possibilities for following the vari- 
ations in solar subsurface features in time with local he- 
lioseismology, provided that the required development of 
efficient and reliable analysis techniques, with the abil- 
ity to keep up with the data flow, takes place. Also, the 
data have great potential for global helioseismology, in- 
cluding the analysis of the near-surface region with high- 
degree modes, once techniques to extract reliable and ac- 
curate frequencies of such modes have been established 
(see Section[3}. 

In the longer term the Solar Orbiter mission is planned to 
observe the Sun from a distance of as little as 0.22 AU 
and, equally importantly, from an orbit inclined at up to 
35° from the ecliptic II71I1 . This will for the first time per- 
mit detailed observations of the solar polar regions, in- 
cluding the possibility of using local helioseismology to 
measure flows at high latitude, such as the near-surface 
meridional circulation which may be crucial for the ac- 
tion of the solar dynamo (Rempel, these proceedings). 
Also, since the spacecraft is orbiting the Sun it will allow 
time-distance seismology with a very long baseline, and 



hence providing information about deep layers in the Sun, 
including the tachocline region. However, the large dis- 
tance from Earth and the resulting low data-transmission 
rate will be a limiting factor for such investigations. It is 
hoped to launch Solar Orbiter in 2015, although this will 
depend on the evolution of the ESA Science Programme 
in a situation of serious financial constraints. 

The DynaMICS plans outlined by Turck-Chieze et al. 
(these proceedings) intend to combine several instru- 
ments for the study of the solar interior and atmosphere 
with the goal to understand, and possibly predict, the 
long-term trends in solar variability. The proposed in- 
struments include the GOLFNG, a further development 
of the GOLF instrument on SOHO; a prototype of the in- 
strument is being constructed and will be tested within 
the next year in Tenerife. The goal is to carry out more 
detailed studies of the deep solar interior, ideally through 
the observation of g modes, and providing information 
about the dynamics and possible magnetic fields of the 
radiative interior. 

The WIRE and MOST missions have demonstrated the 
potential for asteroseismology from space, based on in- 
tensity observations. This potential will be further ex- 
ploited with the launch of the CoRoT satellite towards the 
end of 2006. In particular, the CoRoT long runs, lasting 
5 months, will yield asteroseismic data of unprecedented 
quality for a variety of stars and hence allow much more 
detailed investigations of such stars than previously pos- 
sible. Additional data for a large number of stars will 
be obtained from the Kepler mission, scheduled to be 
launched by NASA in late 2008 id. Kepler will observe 
a field with an area of around 100 square degrees in the 
constellations Cygnus and Lyra continuously for at least 
4 years. The mission is designed to investigate extrasolar 
planetary systems, and in particular search for Earth ana- 
logues, by means of transit detections. To optimize the 
chances of detection of transits more than 100 000 stars 
will be observed at a cadence of 30 minutes. However, 
observations of 512 stars will be carried out at a cadence 
of 1 minute, thus, for example, allowing the detection and 
study of solar-like oscillations. The field contains thou- 
sands of stars for which it is estimated that such oscilla- 
tions can be investigated. The selection of stars can be 
changed every three months; however, it is expected that 
some stars will be observed throughout the mission, to 
improve the quality of the data and search for possible 
frequency variations associated with cyclic stellar activ- 
ity. In addition, the observations with 30-minute cadence 
will obviously be sufficient for the study of solar-like os- 
cillations in giant stars, as well as many other types of 
stellar variability. The Kepler dataset will undoubtedly 
become a unique resource for the study of almost any 
kind of small-amplitude stellar variability. 

Beyond Kepler, there are plans to submit a proposal for 
the PLATO project (PLAnetary Transits and Oscillations 
of stars) for the expected ESA call for new mission pro- 
posals. PLATO will observe a field 10 times as large 
as Kepler, with the goal of detecting and studying large 
numbers of planetary systems as well as carrying out as- 



teroseismic investigations of a very large number of rela- 
tively bright stars. 

Experience from helioseismology has shown the impor- 
tance of obtaining data of the highest possible precision 
and extending over very extended periods. Indeed, it is 
obvious from the analysis reported by Houdek & Gough 
(these proceedings) that a frequency accuracy of the order 
of 0.1 /iHz will be required, for example, for the anal- 
ysis of the signatures of acoustic glitches. Also, it is 
highly desirable to extend the data towards lower fre- 
quencies; here the lifetime of the modes is generally 
longer, leading to more precise frequency determinations, 
the data are more sensitive to the effect of helium ioniza- 
tion, and the frequencies are less affected by the near- 
surface problems. Improving the frequency precision de- 
pends on increasing the signal-to-noise ratio in the ob- 
servations. This evidently requires minimizing the in- 
strumental noise; however, an ultimate limit to the noise 
level is provided by the stellar 'noise', i.e., the effects 
of phenomena in the stellar atmosphere other than the 
oscillations. Already Harvey [57] pointed out that such 
stellar noise is far more important, relative to the oscilla- 
tion amplitudes, in intensity than in velocity observations. 
The results from the VIRGO and GOLF instruments on 
SOHO have clearly demonstrated this, the signal-to-noise 
ratio being higher by almost two orders of magnitude in 
velocity than in intensity. Thus it is highly desirable to 
observe the oscillations in velocity. 

For such observations there is not much to be gained 
from observing from space. However, the obvious chal- 
lenge is to obtain the required continuous observations 
extending over long periods. With existing observa- 
tional facilities this require coordination between differ- 
ent time-allocation committees; also, the required very 
stable spectrographs tend to be on large telescopes where 
obtaining more than a week's observing time is ex- 
tremely difficult. Thus dedicated facilities are needed. 
Nearly continuous observations from a single site can 
be obtained from Antarctica. This is the basis for 
the SIAMOIS project (for Sismometre Interferentiel A 
Mesurer les Oscillations des Interieurs Stellaires) to be 
set up at Dome C, a near-optimal site for astronomical 
observations (Mosser, these proceedings). Here a duty 
cycle of over 90 % is expected during three months in 
mid-winter. 

An alternative is to follow the lead of the helioseis- 
mic networks, such as BiSON and GONG, and estab- 
lish an asteroseismic network of telescopes dedicated to 
Doppler-velocity observations of stellar oscillations. To 
investigate this possibility we have established the SONG 
project (for Stellar Oscillation Network Group). 1 The 
project is currently in an initial phase of conceptual de- 
sign for such a network, including preliminary designs of 
the telescopes and spectrographs. The Doppler-velocity 
observations will use iodine cells as reference; it is esti- 
mated that with one-meter-class telescopes solar-like os- 
cillations can be studied with good signal-to-noise ratio in 
stars brighter than 6th magnitude. Ideally 8 nodes will be 

'see |http : / / astro . phys ■ au ■ dk/SONG| 



required, 4 in each hemisphere, to ensure good coverage 
of the entire sky. The intention is to make the nodes fully 
robotic. The nodes will be established at existing obser- 
vatory sites, to ensure that the required infrastructure is 
available and that limited maintenance can be obtained 
on site. 

It should be noticed that the SONG project and the space- 
based photometric missions are fully complementary. 
Photometric observations can be made simultaneously 
of a large number of stars, providing excellent statistics 
of stellar properties. However, only through Doppler- 
velocity observations dedicated to carefully selected stars 
and carried out over long periods can we achieve the full 
potential of asteroseismology to probe the details of the 
physics of stellar interiors. The combination of these 
techniques, and the parallel development of stellar mod- 
elling, will surely take our knowledge of stellar structure 
and evolution, and of the extreme physical processes that 
control stellar interiors, to new levels of understanding in 
the coming years and decades. 
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